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Abstract 

We describe transverse spectra as well as azimuthal anisotropy ( 112 ) of 
charged hadrons stemming from various centrality PbPb collisions at s/s = 
2.76 ATeV analytically in a ‘soft + hard’ model. In this model, we propose 
that hadron yields produced in heavy-ion collisions are simply the sum of yields 
stemming from jets (hard yields) and yields stemming from the Quark-Gluon 
Plasma (soft yields). The hadron spectra in both types of yields are approx¬ 
imated by the Tsallis distribution. It is found that the anisotropy decreases 
for more central collisions. 


1 Introduction 

Because of its short lifetime, the only way to examine the Quark-Gluon Plasma 
(QGP) formed in ultra-relativistic heavy-ion collisions (HIC), is looking at the par- 
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tides stemming from it. Spectra, angular correlations and their dependence on the 
circumstances of the collision can then be studied. These distributions are effected 
by hadron yields stemming not only from the QGP (we refer to as ‘soft’ yields), but 
also from jets (we call ‘hard’ yields). 

As a first approximation, we make out hadron spectra in HICs as 


ndN n diV hard 
'd^ =P d^ 


dN soft 
d 3 p 


( 1 ) 


where we describe both types of yields by a Tsallis distribution with different pa¬ 
rameters for the following reasons. 

Hard yields: On the experimental side, the Tsallis distribution describes mea¬ 
sured transverse spectra of charged and identified hadrons in proton-proton 
collisions m-m On the theoretical side, the Tsallis distribution provides 
a reasonably good approximation for the transverse spectra of charged pions 
stemming from pp collisions [T3J [14] , and central as well as peripheral PbPb 
collisions obtained via perturbative quantum chromodynamics (pQCD) im¬ 
proved parton model calculations for transverse momenta p t > 4-6 GeV/c 

m- 

Soft yields: The Tsallis distribution has been widely used for the description 
of hadron yields stemming from the QGP |19|-|32|. However, in those models, 
the hard part of the spectrum has not been subtracted. For the emergence of 
the Tsallis distribution in the soft part of the spectrum, there is a chance to 
bring statistical arguments based on non-extensive thermodynamics mm. 
or on super-statistics [7] [221 US [221125] [321 1321 140] . 

We note that transverse spectra and of various identified hadrons measured 
at RHIC energy have been described by a similar model [23j[24l;25j. In that 
model, spectra measured in pp collisions have been used as hard yields, and 
it has been conjectured that hard yields are suppressed at \o\n pr- 

In Sec. [21 analytic formulas are decuced for the hadron spectrum and V 2 - Sec. [3] 
contains fits to charged hadron spectra and V 2 measured in various centrality PbPb 
collisions at yfs = 2.76 ATeV by the CMS |4lll42| and the ALICE [43] collaborations. 
Summary is given in Sec. [4] 


2 Transverse Spectrum and v n 


In statistical models, we obtain the transverse spectrum as a sum of hadrons with 
momentum coming from sources flying with velocities u M as 


dN 

d 3 p 


Too 2tt 


y =o 


= / d< * I da ^ u ^ ■ 


( 2 ) 


-oo 0 
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Here, a is the azimuth angle and ( = | ln[(i + z)/(t— c)]. We parametrize hadron 
momenta as 


p ^ = (tot cosh y , tot sinh y , j»t cos pr sin p) , ( 3 ) 

with y = i ln[(p° +p 2 )/(p° — p z )] and y? being the azimuth angle of the hadron 
momentum. We parametrize the flow as 

u ^ = (7 cosh £, 7 sinh £, 71; cos a, 71; sin a) , ( 4 ) 

with 7 = 1/-v/l — v 2 , and assume that v depends only on a. Though, it is assumed 
that in each source, the momentum distribution of hadrons / is a function of the 
co-moving energy 


UpP /J, \ y=0 ='Y[m T coBh.(;-vpTCOs(<p-a)] , (5) 

the sources may be fireballs [34]—[36], clusters [8j [9 12 37 38 or even jets |39II40| . 
We write the transverse flow as a series, 

OO 

v(a) = vo + 8v m cos (too;) = Vq + Sv(a ) , (6) 

m =1 

and suppose that 8v(a) << 1. We use the Taylor expansion 


/K-p%= 0 


E 

m =0 


[Sv(a)] m 8 m 
1 ! dvH 


f [utff 


v(a) = do 
2 / — 0 


( 7 ) 


and keep only the leading non-vanishing terms in 8v(a). 

Provided that / is a rapidly decreasing function, we approximate integrals with 
respect to ( and p by the maximal value of the integrands times the integration 
interval. Thus, the p integrated transverse spectrum becomes 


dN 

2Trp T dp T dy y = 0 


2t r 



0 


dN 

d 3 p 


y = o 


OO 


E 

m —0 


am d m 

to! dv™ 


f[E(v 0 )\ 


f[E(v 0 )] + 0(8v 2 ) , 


( 8 ) 


with E(vo) = 7o (m T - v 0 Pt ) and a m 

anisotropy becomes 


2tt 

f da [8v(a)] m . Similarly, the azimuthal 

0 


v n = 


2 ’ r dN 

f dpcos(np) p° 

0 


d 3 p 


,=o 8v n jo {vo m T ~ Pr)f [£(> 0 )] 


2 5 „ dN 

2 

f dp p -Ml 
0 dp 

v =0 


f[E{v 0 )] 


0(8v 2 ) (9) 


with Sv n defined in Eq. (JBJ■ 
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For example, in the case of the Boltzmann-distribution / ~ exp [ — /3E(vq)] , 
the anisotropy is 

^BG w Sv nPl 0 (, pT _ VQ TOr ) + 2) (10) 


Thus, oc pt if Pt » m. 

In the case of the Tsallis distribution, / ~ [1 + (g — 1) /3E{yo)]~ 1 ^ q ~ 1 \ the 
anisotropy 


ts _ 3vn f3 _ p T - vp m T _ 

2 1 + (q - l)^7o(mr — vqPt) 


0(6v 2 ) . 


( 11 ) 


Thus, v^ s saturates when (g — 1) /3 70(1 — Vq)pt » 1 - 


3 Fits to the spectrum and V 2 of charged had¬ 
rons in PbPb collisions at snn = 2.76 ATeV 


As conjectured in Sec. [l] we make out the transverse spectrum of charged hadrons 
by the sum of hard and soft yields 


dN 


2 -KpTdpTdy 


y =o 


Ea 


1 + 


(n -1) 

Ti 


[7 i{mT - Vipr) 



-!/(«-!) 


( 12 ) 


(i = soft or hard) where both contributions are assumed to be Tsallis-distributions. 
These yields have maxima at = 7 iUiVi- As long as these maxima are below 
the measurement range, which is the case in this analysis, the isotropic part of 
the transverse flow, v t (denoted by vq in Eq. (JBJ in Sec. [2]) cannot be determined 
accurately. As the dominant part of charged hadrons consists of pions, the argument 
in Eqs. (fl2|l may be approximated by [7,(7717’ — Vtpr ) — m\/Ti « pt/t[ > ° pp with 
the Doppler-shifted parameters 


rjiDopp _ rp 


1 + Vi 
1 — Vi 


( 13 ) 


As can be seen in the top panels of Fig. |T] Eq. (fT2l) describes CMS eh and 
ALICE [43J data on transverse spectra of charged hadrons stemming from PbPb 
collisions of various centralities. Fitted parameters are enlisted in m and shown 
in the bottom panels of Fig. [l] The dependence of the q and T Dopp parameters of 
the soft and hard yields on the event centrality (number of participating nucleons 
Npart ) can be fitted by 


Qi — Q 2 ,z + Pi lB(A^, ar i/2) , 

t d°pp = TTi+rMNpart) ■ (14) 


Though the actual value of the transverse flow velocity cannot be determined in 
this model from the spectra of charged hadrons, it may be guessed using the value 
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Figure 1: Top, transverse spectra of charged hadrons stemming from various 
centrality PbPb collisions at y/s = 2.76 ATeV measured by the CMS [41] (left) 
and ALICE [43] (right) Collaborations. Curves are fits of Eq. dT51> . Bot¬ 
tom, centrality dependence of the fitted q (left) and 5v so ft (right) parameters. 
Straight lines are in Eq. (fl4l) with parameters enlisted in [18]. 


5 










of the QGP-hadronic matter transition temperature obtained from lattice-QCD cal¬ 
culations. As the values of fitted Tffff scatter around 340 MeV, in case of a flow 
velocity of v so ft « 0.6, the real T so f t values would scatter around 170 MeV, which 
is close to the lattice result obtained e.g. in [44] . 

While the tendencies of how fit parameters depend on N part are similar, they are 
not the same within errors in the case of CMS [41] and ALICE [43] measurements. 
It is to be noted that in ra, centrality is determined using the distribution of hits in 
the VZERO detector, which has a rapidity coverage of 2.8 < rj < 5.1 and -3.7 < rj < 
-1.7. In the meanwhile, in ED 02, the collision event centrality is determined from 
the event-by-event total energy deposition in both Hadron Forward calorimeters 
having rapidity coverage of 2.9 < \r]\ < 5.2. 

As seen from Sec. [2] up to O ( Sv 2 ), the transverse spectrum in Eq. (fT2l) results 
in an azimuthal anisotropy of 


V2 


^hard fhard “ 1 “ Wsoft fsoft 
fhard fsoft 


(15) 


where the coefficient functions are 

fivi'y? PT-Vim T 

Wi = -1- ■ ( 16 ) 

1 1 + — [li{rn T - Vip T ) ~ m] 

-L i 

Again, i = soft or hard, Vi are the isotropic part of the transverse flow (denoted by 
v 0 in Eq. ([6]) in Sec. [2]). And Svi are the coefficients of cos(2a) (denoted by Sv 2 in 
Eq. |[6|) in Sec. [2]). 

Fits of Eqs. (fl5ll - (fl6[) to CMS data [42] on V 2 are found in Fig. [2] The 
four different methods used in m for the extraction of V2 are the 2 nd and 4 th 
order cumulant methods denoted by ^{2} and ^{4}, the event-plane V 2 {EP} 
and Lee-Yang zeros V2{LYZ} methods. Fitted parameters are listed in [18] . 

Finally, all four methods for the extraction of V2 in [42] suggest that 5 v so ft 
(the 2 nd Fourier components of the transverse flow of the soft yields) decreases for 
more central collisions (see bottom-right panel of Fig. [I]). This observation is in 
accordance with smaller anisotropy in more central collisions. 


4 Summary 

In this paper, we have simultanously reproduced the transverse spectra and the az¬ 
imuthal anisotropy (i^) of charged hadrons stemming from various centrality PbPb 
collisions at i/s = 2.76 ATeV. In the proposed model, the hadron spectrum is as¬ 
sumed to be simply the sum of yields originated from ‘soft’ and ‘hard’ processes, 
Eq. (jXj) . It is conjectured that hadrons are distributed according to the Tsallis distri¬ 
bution in both types of yields. As for the hard yields, this assumption is supported 
by the observation that the Tsallis distribution provides a reasonably good approxi¬ 
mation for pion spectra obtained via pQCD-improved parton model calculations for 
central or peripheral PbPb collisions at LHC energy [18]. Furthermore, the Tsallis 
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Figure 2: Fits of Eqs. (HSl) (fl6l) to CMS data on vi @2] in the case of four 
types of methods (event plane V 2 {EP}, 2 nd and 4* h order cumulant i> 2{2} and 
i; 2{4} and Lee-Yang zeros V 2 {LYZ} methods). Fit parameters are plotted and 
enlisted in [IS] . 
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distribution describes hadron spectra in pp collisions as well. The soft yields (which 
we identified by what remains of the hadron spectra after the subtraction of the 
hard yields) can also be described by a Tsallis distribution with different parameters. 

Analytic formulas have been obtained for the spectra and for v-i in the limit of 
small transverse flow velocity fluctuations as a function of the azimuth angle. The 
parameters of the soft and hard Tsallis distributions have been determined from fits 
to transverse spectra and V 2 data measured by the CMS [[41] [42] and ALICE [43] 
collaborations. The dependence of the fitted parameters on the event centrality 
( N P art ) have been found similar in the case of the CMS and ALICE data. Fits 
to CMS data on V 2 suggest that in this model, the anisotropy decreases for more 
central collisions. Fit parameters are enlisted in [Si- 
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